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ABSTRACT 

Doppler spectroscopy of p 1 Cnc has detected evidence of a companion with an orbital period 
of 14.65 days and a minimum mass of 0.88 Jupiter masses. Astrometric observations performed 
with the Hubble Space Telescope Fine Guidance Sensor lr using a novel new observing technique 
have placed an upper limit on the astrometric reflex motion of p 1 Cnc in a time period of only 
one month. These observations detected no reflex motion induced by the 14.65 day period radial 
velocity companion, allowing us to place a 3cr upper limit of ~0.3milliarcseconds on the semi¬ 
major axis of this motion, ruling out the preliminary Hipparcos value of 1.15milliarcseconds. 
The corresponding upper limit on the true mass of the companion is ~30 Mj, confirming that it 
is a sub-stellar object. 

Subject headings: astrometry—planetary systems—stars: low-mass, brown dwarfs 


1. Introduction 

The Doppler spectroscopy technique has been 
used in recent years to detect low amplitude, peri¬ 
odic radial velocity variations in ~60 nearby stars, 
which have been interpreted as due to planetary 
mass companions. This interpretation requires 
several assumptions, namely that the root cause 
of the variation is Keplerian in nature; that the 
companion mass (M c ) is substantially less than the 
mass of the primary; and that one is seeing light 
from a single star as opposed to an unresolved, 
comparable-mass binary system (Imbert & Prevot 
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1998). The quantity determined by the radial ve¬ 
locity observations is then M c sin(i) where i is the 
unknown inclination of the orbital plane to our 
line of sight to the star. The argument usually ad¬ 
vanced to suggest that the masses of the compan¬ 
ions must be small assumes that the distribution of 
orbital inclinations is uniform so that, on average, 
(M c ) = 4(M C sin(i))/7r (Chandrasekhar & Munch 
1950), which implies that the true companion mass 
cannot be much greater than the minimum mass, 
M c sin(i). Note that this statement is only true 
for a large sample, and it does not preclude indi¬ 
vidual companions from having true masses that 
are significantly larger than their minimum mass. 

Several pieces of work support the planetary 
mass interpretation. For example, Ito & Miyana 
(2001) have determined that the dynamical sta¬ 
bility of the Ups And system constrains the mass 
of the outer planet to be less than 1.43 times 
its minimum M c sin(?) value. There is also one 
case, HD 209458, for which the companion has 
been observed to transit the disk of the central 
star (Charbonneau et al. 2000), hence giving a 
measure of both the companion radius and the 
inclination of its orbit. The companion mass of 
~0.6 Jupiter masses (Mj) is consistent with that 
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of a giant planet. However, despite the statistical 
and dynamical arguments, and the special case of 
HD 209458, doubts have persisted that M c siu(i) is 
nearly equal to M c for several reasons. 

Using a combination of Hipparcos and ground- 
based MAP (Multichannel Astrometric Photome¬ 
ter) astrometry in conjunction with the orbital 
parameters derived from the radial velocity (RV) 
data, Gatewood, Han & Black (2001) have found 
that the M c sin(i) = 1 Mj companion to p CrB is in 
a nearly face-on orbit with i = 0.5°, and has a true 
mass of 115 Mj, making it an M dwarf rather than 
a planetary-mass object. Zucker & Mazeli (2000) 
have used Hipparcos astrometry to show that the 
M c sin(i) = 6.4 Mj companion to HD 10697 is ac¬ 
tually a 38 Mj brown dwarf in an orbit with an in¬ 
clination of just 5°. A similar analysis combining 
the Hipparcos and radial velocity data extended to 
the 30 systems with orbital periods in excess of ten 
days (Han, Black & Gatewood 2001) suggests that 
at least four of the 30 stars they analyzed have stel¬ 
lar mass companions, that is, M c > 80 Mj, which 
in turn would require i < 1°. If the distribution 
of inclinations in the sample is uniform, the prob¬ 
ability of a system having an inclination i < i a is 
given by 1 —cos(i 0 ). For i a = 1°, the probability 
would be 1.5 x 10 -4 , making it unlikely that even 
one such system would be observed in a sample of 
2000 stars. This led Han et al. to suggest that 
there might be a bias toward small inclination an¬ 
gles in the radial velocity studies. Pourbaix (2001) 
also finds statistically significant astrometric or¬ 
bits with low inclinations for three out of four of 
the Han et al. potential stellar mass companions. 
However, Pourbaix & Arenou (2001) argue that 
the trend to low inclinations is an artifact of the 
adopted reduction procedure, and that the astro- 
metric data are not precise enough to allow the 
conclusion that a significant fraction of the RV 
companions have stellar masses. 

An independent line of evidence also raises the 
possibility of stellar mass companions for some of 
the 60 systems. Among a subset of 9 of the stars 
of spectral type F with candidate planetary com¬ 
panions, Suchkov & Schultz (2001) have identified 
three potential binaries, HD 19994, HD 89744, and 
HD 169830, along with HD 114762 (private com¬ 
munication) , based on the fact that they are overly 
bright for their spectral type and Hipparcos dis¬ 
tance. Gonzalez et al. (2001) found that two stars 


with RV companions, HD 37124 and HD 46375, are 
similarly too luminous for their spectral type and 
distance; they might each be unresolved binary 
systems. However, they also note that a wide, 
long-period binary cannot be ruled out, so the 
short period RV companion need not be the source 
of the “excess” luminosity. 

Another line of reasoning suggesting that the 
RV companions may not be planetary in nature 
comes from an analysis of the distribution of their 
eccentricities and orbital periods, which are statis¬ 
tically indistinguishable from those for single-line 
spectroscopic binaries (SBls) (Stepinski & Black 
2001; Heacox 1999). Moreover, the apparent cor¬ 
relation of eccentricity with orbital period for the 
RV companions is similar to that for SBls (Black 
1997; Heacox 1999), and their bivariate probabil¬ 
ity distribution functions are again statistically in¬ 
distinguishable (Stepinski & Black 2001). The ob¬ 
served orbital properties also differ strongly from 
those of our own planetary system. The current 
data, interpreted as planetary systems with ran¬ 
dom inclinations, also does not lead to a simple 
theoretical picture (e.g., Marcy et al. 1999), and 
there is currently no compelling dynamical argu¬ 
ment to support the interpretation that all of the 
RV companions are planets, rather than brown 
dwarfs or stars in low inclination orbits. A de¬ 
termination of which companions (if any) are not 
planetary would facilitate the development of a 
dynamical model. 

Astrometric observations of the stellar reflex 
motion induced by a companion can potentially re¬ 
move the uncertainty in sin(i), and thus determine 
the companion mass. The Hubble Space Tele¬ 
scope (HST) Fine Guidance Sensor lr (FGSlr) 
can measure relative stellar positions to an accu¬ 
racy ~0.3milliarcsecond (hereafter mas), a factor 
of 3-5 improvement over that of Hipparcos and 
MAP data. In this paper we present the results 
of a pilot program utilizing an observing tech¬ 
nique designed to quickly search for perturbations 
larger than about 0.3 mas to the position of p 1 Cnc 
(HR3522, HD 75732, 55 Cnc), a G8V star with an 
M c sin(i) = 0.88 Mj companion in an orbit with 
a period of 14.65 days (Butler et al. 1997). Han 
et al. report a preliminary reflex motion with a 
semi-major axis of 1.15 mas for this star. We dis¬ 
cuss our rationale for selecting this target and our 
observing strategy in section 2, and our results in 
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section 3. 

2. Observations 

Detection of companion-induced reflex motion 
is complicated by the need to determine the typi¬ 
cally much larger proper motion (p a , ps) and par¬ 
allax ( 7 r) of the star, which normally requires ob¬ 
servations over a baseline of at least one year. 
However, a carefully designed observing strategy 
can shorten this time dramatically provided an 
RV system with a relatively short period is cho¬ 
sen, and provided the individual astrometric ob¬ 
servations can be made with sufficient precision. 
We therefore concentrated on stars from Han et 
al. with the following characteristics: relatively 
short period; a large (> lrnas) reflex motion in¬ 
ferred to be present with at least moderate sta¬ 
tistical significance; and ones most likely to have 
companions that are M dwarfs or brown dwarfs 
(their Groups 2 and 3). An additional criterion 
was the suitability of the target for FGSlr astrom¬ 
etry, namely, the availability and distribution of 
stars in the FGS field of view, needed to define 
the inertial reference frame. Using these criteria, 
p 1 Cnc was determined to be the most suitable tar¬ 
get. Its predicted reflex motion of 1.15 mas ra¬ 
dius implies that the companion is an M-dwarf of 
126 Mj rather than the M c sin(i) value of 0.88 Mj. 
However, the Han et al. result has only borderline 
statistical significance for this target, and there¬ 
fore falls into their Group 2, stars for which they 
expect the majority of the companions to be brown 
dwarfs. 

We realized that we would detect a reflex mo¬ 
tion with FGSlr only if the RV companion is more 
massive than about 40 Mj which, to be consistent 
with the radial velocity data, implies that the in¬ 
clination of its orbit would be nearly face on. (In¬ 
clinations larger than about 2° imply a companion 
mass that is too small to produce an astrometric 
signature large enough to measure with HST.) Our 
observing strategy was designed to optimize our 
ability to detect a perturbation with the known 
14.65 day period and low eccentricity of the com¬ 
panion’s orbit, and the unknown value of H (the 
longitude of the ascending node). 

We performed a set of FGSlr observations in 
the one month period centered around the time of 
maximum parallax factor in right ascension (RA) 


(hereafter a max ) on 1 May 2001. We observed at 
pairs of epochs with times that were both phased 
with the companion’s 14.65 day period (P), and 
symmetric about the time of a max . Explicitly, ob¬ 
servations were executed in pairs occurring at the 
same phase from April 17 through May 16. The 
relative phases for each epoch are shown in Fig¬ 
ure 1. With additional observations on May 30 
(epoch 10 at +2P), we monitored the star’s po¬ 
sition over more than two full orbits of the com¬ 
panion. The star’s proper motion p a can then be 
determined from the epoch pairs without a simul¬ 
taneous determination of its parallax. Since the 
rate of change of the parallax as projected along 
declination, dns/dt, was constant to within 1% at 
the time of a raffl , it simply added to p$ as a con¬ 
stant. Measurable reflex motion would manifest 
itself as a larger dispersion in the proper motions 
of the pairs. This strategy favors shorter period 
systems since the reflex motion is then more ap¬ 
parent, i.e., less diluted by the star’s much larger 
proper motion. 

Epochs 2, 4, 5, 6, and 8 each consisted of two 
HST visits, the other epochs of 1 visit. The data 
obtained in epoch 3 were severely degraded due 
to a loss of lock on guide stars part way through 
the observing sequence, so the proper motion mea¬ 
surement utilizing data from the epochs 3 (— 4P) 
and 7 (+(jP) temporal pair could not be made. 
Each visit consisted of one HST orbit. All obser¬ 
vations were obtained at a fixed HST roll angle 
and pointing, and the observing sequence for each 
visit consisted of 4-5 observations of p 1 Cnc inter¬ 
spersed with 1-4 observations of each of the 7 ref¬ 
erence stars. One of the reference stars, pCncB, 
is the proper motion companion to p 1 Cnc. Obser¬ 
vations of this star were unsuccessful in the first 6 
epochs due to incorrect input coordinates. 

3. Results 

In a given HST orbit the relative positions 
of the stars were determined with a precision of 
about 1 mas, as illustrated by the residuals for 
p 1 Cnc shown in Figure 2. Data from multiple 
epochs were combined by a standard 6 parame¬ 
ter plate solution using the GAUSSFIT program 
(Jeffereys, Fitzpatrick, & McArthur 1987). The 
solution yields the proper motion of p 1 Cnc and 
p 1 CncB, and a catalog of star positions in (£,?y) 
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space, along with residuals to these quantities. 
When combining data from visits restricted to the 
± temporal pairs defined above, the parallactic 
displacement along RA is not an issue since it is 
common to such data sets. However, when other 
visits are included in the analysis, the parallax 
of p 1 Cnc must be accommodated in the solution. 
Since it was not possible to solve for the paral¬ 
lax (because of the short interval over which the 
observations were made), we adopted the Hippar- 
cos value of 7 r = 79.8 mas for p 1 Cnc and pCncB, 
while constraining the parallax and proper motion 
of the remaining reference stars to zero. To deter¬ 
mine our sensitivity to the choice of n, solutions 
were done with 75 < n < 82 mas. The residu¬ 
als in the astrometric catalog produced for each 
assumed value of i r were essentially constant, in¬ 
dicating that our particular choice of 7 r was not 
important provided it was correct to within a few 
percent. 

The residuals of the star positions in the derived 
astrometric catalog were small, 0.27-0.45mas 
for £, and 0.28-0.60 mas for 77 , (see Figure 2) 
with the exception of those for pCncB, which 
were ~1 mas, attributable to the fact that it 
was observed in only 5 of 15 HST orbits. The 
derived proper motions for p 1 Cnc and p Cnc B 
were ( p a ,ps) = (—515.2 ± 6.6,-256.6 ± 4.7) 

and (-508.7 ± 18,-221.7 ± 12.6) mas/yr re¬ 
spectively. The error bars include the cumula¬ 
tive 1 -sigma errors from all the fitted parame¬ 
ters. This compares with the Hipparcos values 
of -485.46 ± 1.03, -234.40 ± 0.72 for p 1 Cnc. The 
difference between our proper motion and the Hip¬ 
parcos value is most likely due to unmodeled, but 
small, proper motion of one or more reference 
stars, which can’t be accurately measured over 
such a short time. We emphasize that our tech¬ 
nique and results do not depend on an accurate 
determination of the absolute proper motion, but 
rather on the differential accuracy with which we 
can determine it among epochs and epoch pairs, 
so a systematic offset from the Hipparcos value is 
not relevant. 

We used three different techniques to search 
for reflex motion and to gauge the sensitivity of 
our measurements. First, we fixed the inclina¬ 
tion at i = 0 ° and the period of the companion 
to P=14.65days and solved for the phase and ra¬ 
dius of the reflex motion, along with residuals of 


these values, using data from all 10 epochs. This 
solution yielded a reflex motion with a radius of 
0.08 ± 0.2 mas, i.e., a non-detection. 

Second, we again fixed the inclination and pe¬ 
riod of the companion, then in simulations im¬ 
posed a range of values for the phase angle and 
radius of the reflex motion on the data, and solved 
for the phase and radius (again using data from 
all 10 epochs) to assess how accurately we could 
recover them if such perturbations were present. 
The recovered values of {r,<f>) matched the imposed 
values to within the residuals computed by the 
model for r ;> 0.3 mas. From this test we con¬ 
cluded that our observations would have detected 
any 14.65 day period circular reflex motion with 
a radius in excess of 0.3 mas at 1.5 a. No such 
motion was seen in the real data. 

Our final and most sensitive test for detect¬ 
ing the reflex is to look for dispersions in the 
proper motions computed from the ± tempo¬ 
ral pairs taken at epochs with identical paral¬ 
lax factor a. The proper motion computed from 
epochs (1,9) and (3,7), each pair taken at the 
same star/companion orbital phase, should agree 
and be the real proper motion. Unfortunately, 
as mentioned earlier, the astrometry obtained in 
epoch 3 were compromised, so only the data from 
epochs (1,9) could be used to measure a “reflex- 
free” proper motion. Table 1 (“Measured”) shows 
the values of (p a ,ps) computed from the three 
remaining ± temporal pairs. As expected, the av¬ 
erage value of (p a ,ps)—( 516,251) mas/yr derived 
from the pairs is very close to that derived using 
the data from all 10 epochs. However, the disper¬ 
sion of ps is nearly a factor of three larger than 
that of p a , due primarily to the (1,9) pair, which 
we attribute in part to the fact that the epoch 1 
and 9 observations consisted of only a single HST 
orbit, while the epoch 2, 4, 6 , and 8 observations 
consisted of two HST oribts each. 

Note that the change of the star’s position due 
to its orbit around the star-companion barycen- 
ter has a position angle which is 180° different in 
epochs (2,8) relative to that in epochs (4,6) (see 
Figure 1). Therefore any astrometric signature of 
the star’s companion would manifest itself most 
dramatically in the the proper motions computed 
from these epoch pairs. The data from epochs 
(4,6) take on particular significance because they 
span the shortest time interval among the tem- 
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poral pairs, hence any detectable reflex motion is 
least diluted by the star’s true proper motion com¬ 
puted from this pair. 

To illustrate this, we conducted a series of simu¬ 
lations whereby we again impressed upon the data 
a circular reflex motion with a 14.65 day period, 
along with semi-major axes of 0.3, 0.5 and 1.0 mas, 
and a variety of phase angles. The results of these 
simulations are shown in Table 1. As expected, 
because it is the “reflex free” pair, the epochs 
(1,9) proper motion shows no change due to the 
imposed perturbation, while the (4,6) pair has 
changed significantly, even for the smallest 0.3 mas 
perturbation. In these simulations, a reflex mo¬ 
tion is absorbed into the proper motion computed 
from each ± temporal pair. Therefore, the greater 
the reflex, the larger the dispersion in the val¬ 
ues of (p ai ps), even though the perturbation has 
little effect on the average value of the proper 
motion computed from all the temporal pairs. 
For example, the measured dispersion in (p a ,p$) 
is (5,14), while that from the simulation with 
(r, 4>) = (0.3 mas,103°) is (18,15) and that from 
(r, (f>) = (0.3 mas,13°) is (5,25). We take the ratio 
between the standard deviation of the measured 
proper motion and the perturbed proper motion 
as a rough (and very conservative) estimate of our 
ability to detect a real perturbation to the proper 
motion. The standard deviations in Table 1 are 
computed in the conventional fashion, i.e., as the 
square root of (valuer — average) 2 divided 

by (n — 1) using the three epoch pair values listed 
in the Table for the three cases presented. For 
the 0.3, 0.5 and 1.0 mas perturbations these ratios 
(average for the two phases presented in Table 1) 
are 3, 4, and 7. Comparison of the measurements 
and simulations for only the (2,8) and (4,6) pairs 
would in fact be a more valid (and less conserva¬ 
tive) measure of the sensitivity because the (1,9) 
pair is unaffected by the perturbation, and would 
give an even larger ratio. We therefore conclude 
that a reflex motion of 0.3 mas is conservatively 
ruled out at about the 3<r level, and any reflex mo¬ 
tion with a semi-major axis in excess of 0.5 mas is 
firmly ruled out at greater than 4 a. 

In summary, the strategy we chose for attempt¬ 
ing to astrometrically detect and measure a reflex 
motion of p 1 Cnc with HST/FGSlr yielded the an¬ 
ticipated sensitivity. With just 14 orbits of HST 
time over a 30 day interval we were able to con¬ 


clusively rule out the preliminary 1.15 mas pertur¬ 
bation proposed by Han et al. Furthermore, we 
ruled out at the 3cr level a perturbation with an 
amplitude greater than about 0.3 mas. This places 
an upper limit of ~30Mj on the mass of p 1 One’s 
companion, implying that it is a sub-stellar ob¬ 
ject. We emphasize that this observing technique 
is fully capable of detecting the larger and more 
statistically significant astrometric perturbations 
for, e.g., pCrB and HD 195019 found by Han et 
al. (2001) and Pourbaix (2001) if they are real. 
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Table 1. Proper Motion Measurements and Simulations 


Measurements [n a , ns 

in units of 

mas/yr] 


Epochs (1,9) 
(2,8) 

. (4 .’ 6) 

Standard Deviation 

-520,-274 

-509,-241 

-520,-239 

5,14 



Simulation 1. Phase= 

103°(RA) 



Radius 

0.3 

0.5 

1.0 

Epochs (1,9) 

-520,-274 

-520,-274 

-519,-274 

(2,8) 

-501,-241 

-497,-241 

-485,-241 

(4,6) 

-551,-233 

-571,-232 

-621,-231 

Standard Deviation 

18,15 

27,15 

50,16 

Simulation 2. Phase= 

13°(Dec) 



Radius 

0.3 

0.5 

1.0 

Epochs (1,9) 

-520,-274 

-520,-274 

-520,-274 

(2,8) 

-508,-248 

-508,-253 

-508,-265 

(4,6) 

-520,-203 

-519,-183 

-518,-132 

Standard Deviation 

5,25 

5,34 

5,56 



-P -3/4P -1/2P -1/4P 0 1/4P 1/2P 3/4P +P 

Time 


Fig. 1.— The relative phase of observations at 
each epoch is illustrated schematically on this plot 
of parallax factor in RA vs time (in units of the RV 
companion period, P). Because we are searching 
specifically for a companion with a circular, face- 
on orbit with a period of 14.65 days, the absolute 
phase of the observations is not important. 
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Fig. 2.— (a) The X and Y residuals for all ob¬ 
servations of p l Cnc relative to the best-fit 6 pa¬ 
rameter plate solution model with proper motion 
and parallax removed. Data from the 10 epochs, 
excluding epoch 3, are plotted using symbols as in 
(b). (b) The £, rj position of p 1 Cnc at each epoch 
from the best-fit model, illustrating the accuracy 
of ~0.3mas achieved by the observations. 
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